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Getting (3 — a without penguins 
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Oscillation effects in B° — > KsD° and related processes are considered to determine S = 
/3 — a + Tv — 2/3 + 7- We suggest that D° decays to CP eigenstates used in concert with inclusive 
D° decays provide a powerful method for determining 8 cleanly i.e. without any complication 
from penguin processes. The CP asymmetry is expected to be<^40% for D° decays to non-CP 
eigenstates and <^ 80% for decays to CP eigenstates. This method can lead to a fairly accurate 
determination of S with O(10 8 - 10 9 ) B -mesons. 
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The two asymmetric B-factories have made remark- 
able progress in determining one of the angles {j3) of 
the unitarity triangle; the world average now stands at 
sin2f3 WA = 0.78 ±0.08 |, §. This is in very good agree- 
ment with the expectations from the Standard Model 
(SM), sin2(3 SM = 0.70 ±0.10 @. However, considerable 
amount of theoretical input has to be used to deduce 
sin2/3 SM and progress in reducing the theory error is 
likely to be rather slow. Thus, methods that determine 
the angles without the uncertainties of hadronic matrix 
elements are crucial in testing the CKM paradigm Q to 
an increasing degree of accuracy in an effort to search for 
CP-odd phase(s) due to physics beyond the SM. 

In the SM, CP violation is controlled by only one CP- 
odd phase. Therefore, different decays which measure 
the same angle of the unitarity triangle (UT) may give 
inconsistent results if physics beyond the SM is present. 
Likewise other apparent failure of unitarity of the CKM 
matrix, such as the failure of the UT to close, would also 
indicate new physics. Beyond the phase (3, the determi- 
nations of a and 7, therefore, provide key SM tests. 

Two extensively studied methods for determining a 
already exist, via B — > 7T7t || and B — > pn [||. In 
these approaches, in addition to some experimental diffi- 
culties, considerable theoretical input is essential as the- 
ses modes receive large QCD-penguin as well as some 
electro-weak penguin (EWP) contributons. While efforts 
at these methods should certainly continue, it is also very 
important that, in our drive towards precision, we de- 
velop methods that require no theoretcal assumption and 
therefore have zero theory error. The key point is that 
effect of beyond the SM CP-odd phase(s) on B-physics 
may be quite small so any residual theory error on the 



determined unitarity angles may mask the effect of new 
physics and thwart experimental searches. 

In this work we wish to report on our study of a method 
to extract 8 = (3 — a + tt = 2[3 + -f that uses interference 
between b — > u and b —> c tree graph exchanges only; 
no penguin contribution, strong or EW, or any theoreti- 
cal assumption is involved. Given that (3 is already well 
measured, this method is very effective in determining 
a "cleanly", i.e. without QCD complications. In ad- 
dition, this method can also be used to simultaneously 
extract (3, allowing a crucial check against the value of 
(3 determined with the B — * J/ipK® approach |l|, ||, 
As mentioned before, a difference in the two determina- 
tions of (3 may then be an indication of new physics. The 
basic idea behind the method has already received some 
attention §, g [fij ~ 
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|l2fl . We extend and complement 
these earlier studies in several ways so that it becomes 
now a powerful approach to determine a, and possibly [3, 
without any complication from penguins. 

In principle, a comparison of time dependent CP asym- 
metry measurement in B°(B°) — ► KsD° with that in 
KsD° suffices to give 6 || |ll| . In practice, though, as has 
already been noted previously, flavor tagging of D°(D°) 
appears extremely difficult jl0|, [lj|. Semi-leptonic tags 
suffer from very serious background from prompt B- 
decays, B — > IvXc] therefore, here we will not consider 
the possibility of semi-lcptonic tags further. Hadronic 
tags of D° (say via D° — » K ~ir + ) receive appreciable cor- 
rections from doubly-Cabibbo-suppressed decays of D°. 
As in the case of 7 extraction with B ± p3| , this inter- 
ference can be used to our advantage in determining S as 
Kayser and London (KL) have discussed p0| . 

In this letter, we would like to highlight at least two 
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additional methods which will be shown to have great 
practical importance in extracting S. First of all, S may 
be determined if D° decays to CP eigenstates (CPES) 
are observed, provided both CP=+1 and CP=— 1 states 
are used. Although we find that while neither this CPES 
method, nor the CP non-eigenstate (CPNES) method of 
KL [^ot , can separately provide an especially sharp deter- 
mination of S, a great improvement is achieved if both 
approaches are used together because both data sets de- 
pend on a common set of parameters in the amplitude. 
Secondly, we will generalize these methods from single fi- 
nal states to inclusive sets of final states. In this way we 
can use the entire observable hadronic branching ratio 
of the D° greatly enhancing the statistical power. Fi- 
nally we will briefly discuss methods whereby ancillary 
information constraining 5 may also be obtained. The 
methods which we describe share with KL|l(| the fea- 
ture that the amplitude parameters are overdetermined 
and therefore a value of j3, in addition to 5, may also 
be extracted from the same data providing a valuable 
comparison to (3 obtained from B — > J/ifrKg. 

Consider now the case where B°(t)/B°{t) -> K S D°, 
K S D° followed by the decay D°/D° -> F; F denotes an 
inclusive set of states F = {fi} and in general F ^ F. For 
example, the set {fi} may range over states of different 
particle content (e.g. K~ + mr) or different points in 
phase space jlj] (e.g. each /, is a point on the K~tt + tt° 
Dalitz plot) or a combination of both. For each fi the 
four relevant amplitudes are: 
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where, we have adopted the Wolfenstein]15[ represen- 
tation of the CKM matrix, and without loss of gen- 
erality, we can choose the strong phase convention so 
that Ai — A is real. The quantity is the ra- 
tio |„4(_D° — > fi)/A(D° —* fi)\ which we will assume is 
known from the study of D° decay. The strong phase 
VD(fi) = arg (A(D° ft/AiD -> /;)) we will assume 
to be not known apriori. Likewise the parameter tb and 
the strong phase t\b given by rge" 12 = e~ n A(B° — > 
K S D°)/A(B° -> K S D°) are also assumed to be not 
known apriori. Note that {r£>, r)u, A} depend on the 
state fi while {tb, t/b} are independent. 

The time dependent decay rates for this decay is: 
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Mfi) + Mfi) and Af 
efficients X, Y and Z in Eqn 

Ei(\A(fi)\ 2 +\A(f, 

and Z(F) = Y, i Im(e-' Mt) A{fi)*A{h)). We can expand 
these quantities in terms of eqn. ([!]) and obtain 
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where r = Tst and xb = Arns/rs while b = +1 
for B(t) and b = -1 for B(t). Defining A(fi) = 



where A 2 = J2 t A 2 (f l ) 1 f 2 D = A 2 (f, 

R F e">° = (EiMfi)rD(fi)e^°^)/{Afg). 

The corresponding quantities for F are given by 
X(F)( V b, Vd, 7) = X(F)(-r) B , -VD,_l); Y(F)( Vb , 
tlD,l) = ~ y ( f )("Vb, -Vd, 7) and z ( f )(vb, Vd,j) = 
Z(F)(—r)B, —?lDi 7) assuming that there is no additional 
CP violation in D° decay filf. 

Initially we will assume that (3 is well determined. Let 
us now consider the special case where F consists of 
CPES with eigenvalue a = ±1. In this case, the modes 
add coherently and so Rf = 1, ?£> = 1 and f/F> = or 
7r for a — +1, —1 respectively. The three observables 
X(F), Y(F) and Z(F) thus depend on the four param- 
eters {A, tb, Vs, 8}. If we have the two data sets, for 
a = +1 and for a = — 1, then there are five independent 
observables (note that Y(a = +1) + Y{a = -1) = 0) 
determining the same four parameters and so the system 
is overdetermined and one may solve for i5. 

Some examples of CP=— 1 final states jlTj include 
K s ir° (BR=1%), Ksv (0.35%), K s p° (0.6%), K s lo 
(1.1%), K s r]' (0.9%) and K S (j) (0.4%) giving a total of 
about 4.4%. CP=+1 final states include K s f (0.3%), 
7T+7T- (.07%) and K+K- (.21%) for a total of 0.6%. For 
each of the modes with a K$ one can construct a mode 
of the opposite CP by changing the Ks to a Kl. We can 
also change the Ks which arises from the B° decay i.e. in 
B° —> KsD° (which we refer to as the fast kaon) to a K^. 
Switching the fast kaon to K~l changes 773 — > t\b + K and 
thus gives the same information as switching the slow 
kaon (i.e. the kaon arising from D° decay). It should 
be emphasized that, in this instance, including the final 
states both with Ks and with Kl, increases the num- 
ber of observables and, as mentioned above, enables the 
system of equations to become soluble. This is in con- 
trast, for example, with the case of B° — * J/ipKs versus 
B° — » J /ipKi, where switching the kaon merely improves 
statistics but does not provide additional independent 
observables. 

We can extend this CPES method to consider inclusive 
final states. If F is defined in a CP invariant manner (eg. 
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F = Ks + X, BR=21%) the resultant observables will be 
similar to the pure eigenstate case. Here again f d = 1 
and f]£> — or tt but Rp , which measures the purity of 
F, will not be 1. The 3 observables are thus dependent 
on 5 parameters {A, rs, f/s, 5, Rf}. As before, we can 
obtain a solution by changing the fast K° to a Kl and/or 
changing the slow K°, in the case where Kg G F. Again, 
these K° changes will lead us to 5 observables. 

In jl(J KL studied the special case where F consists 
of a single quantum state which is a CPNES (e.g. / = 
K~n + ). Then Rp = 1 but {^4, rs, t]b, S, t]o} are not 
known. If we take the point of view that and all the 
relevant D° branching ratios are well determined then 
as discussed in |l(J there are 6 observables {X(f), Y(f), 
Z(f), X(f), Y(f), Z(f)} determining these 5 parameters 
and so the system is overdetcrmincd; therefore, one can 
extract 6. Furthermore, as in |To| one can also take the 
point of view that f3 is a free parameter and solve for 
both P and 5 from the same six observables. In this 
context, as in the CPES case, taking the fast kaon to 
be Kjj (rather than Kg) provides 6 more independent 
observables dependent on the same parameters rendering 
the system even more overdetermined. 

There is a great advantage to combining the CPES and 
CPNES methods above since the parameters involved in 
the CPES case are a subset of those for a CPNES. Thus, 
combining information from CPES and CPNES meth- 
ods can increase the number of observables to nine or 
eleven depending on whether one or both CP eigenval- 
ues are included, respectively. Indeed, if also the fast Kl 
is taken with the CPNES then the number of observ- 
ables increases to seventeen. The number of parameters, 
of course, stays the same, i.e. five (or six if we also in- 
clude f3 as an unknown). Thus, not only there is enough 
information but in fact there is considerable degree of 
redundancy to solve for the unknown parameters. 

Likewise considering several CPNES can enhance the 
degree of over determination. For each CPNES we add, 
we have six new observables but introduce only one new 
parameter (j]d(F)) giving a net gain of 5. Indeed there 
are several candidate modes: K~ir + (branching ratio 
3.8%), K-p+ (10.8%), K*-n+ (5.0%), K*°ir° (3.1%), 
K*-p+ (6.1%) and K*- a + (7.3%) giving a total 36%. 
In this method one would have to separate the quasi two 
body modes from the broad resonances (eg K~ p + ) mak- 
ing it somewhat difficult. 

Generalizing the CPNES case to inclusive states should 
provide the most statistically powerful data to determine 
8. For instance, the inclusive D a — > K ~+X has a branch- 
ing ratio of 53% fT?]] . In this case, we have the general 
case of eqn. (||) and so six observables {X(F), Y(F), 
Z(F), X(F), Y(F), Z(F)} are determined by the six pa- 
rameters {A, rg, t]b, ro, Vd, 8} and so the system can 
be solved with some discrete (8-fold) ambiguities. 

This may be improved in two ways. Firstly, one can 
segregate the set F into several subsets. Thus each 



additional set F provides six new observables but in- 
troduces only two new parameters (rfo and Rf) giv- 
ing a net gain of 4. For instance, a substantial frac- 
tion of K~ + X is made up of the exclusive state K~ir + 
(4%) together with the inclusive (in the sense that these 
modes depend on phase space variables) states K~ir + ir Q 
(13.9%), K-iT+ir+ir- (7.5%), K-ir+n-ir (10.0%) and 
K ~7r + 7r + 7T~7r° (4.0%) giving a total branching fraction 
of about 40%. Another approach would be to divide the 
K~ +X into separate bins according to the energy of the 
K~ in the D° frame. This would approximate the above 
since a higher energy K~ would tend to be associated 
with fewer pions; then one would not have to identify 
the content of the X state. Secondly, one could combine 
the F (inclusive) method with the CPES method. 

The magnitude of the time dependent CP asymmetry 
for various final states can be seen in the expression for 
Z in Eqn. (|). If D° -> F is Cabibbo allowed then r D w 
sm 2 9 c « 0.05 while r B « \V ub \\V cs \/ {\V cb \\V us \) « 0.36. 
The dominant term in Z is thus the fourth term which 
will lead to CP violation of < 36%. Note that this term 
does not become small in the limit Rp — » 0. In the case 
where F is a CP eigenstate, so rp, = 1, the second term 
oc sin2/J « 0.8 becomes dominant. If Rp is small then 
the third and fourth terms are dominant, giving again a 
contribution oc r B « 0.36. 

Now, in order to illustrate the relative power of 
the different methods, let us consider the following 
toy model. First, we estimate BR(B° — > D°Kg) ~ 
sm 2 9 c (l/N 2 )BR(B° -> D~-k + )/2 m 10~ 5 . For this ex- 
ample, we will arbitrarily take rj B = 50° and r\r> = 70° 
with 7 = 60° and (3 — 25°, consistent with the B factory 
values and so S = 110°. 

In Fig. 1 we plot the x 2 which would be obtained as- 
suming that Nb —(number of B mesons) (acceptance)= 
10 9 for various combinations of the data. The thin solid 
line gives the minimum value of % 2 as a function of 5 
obtained for the data from the single final state K~ir + . 
Clearly discrete ambiguities in the solution tend to con- 
spire to keep the value of \ 2 relatively low. The dashed 
line gives the minimum x 2 using CP eigenstates contain- 
ing a Kg and a K^. The dotted line shows the case 
when K~ir + data is combined with CPES containing Kg. 
The dashed-dotted line gives the result for the inclusive 
K ~ + X alone where we have taken for the purposes of 
illustration t)d — 70°, fp> the same as for the K~ir + and 
Rf = 0.1 and the thick solid line combines this with the 
CP= — 1 eigenstates. In all cases, we have assumed that 
the overall tagging efficiency for the B° flavor is 25%. 

Table 1 shows the one sigma error on 6 for various in- 
puts. Clearly, the best results are obtained when the ob- 
servables overdetermine the parameters and a large frac- 
tion of the BR is included in the sample. Thus when 
K~ + X is used together with Kg CPES, error on S is 
±2.5° with N B — 10 9 ; with Nb = 10 s this error increases 
to ±11.4°, so is still quite useful. 
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through £?° — > D°Ks can be greatly enhanced by con- 
sidering Z?° decays to CP eigenstates and by using in- 
clusive sets of D° decays. In particular, using inclusive 
D° decays such as K~ + X together with CP eigenstates, 
our illustrative calculation suggests that as the number of 
available B mesons increases from 10 8 to 10 9 a determina- 
tion of 6 with a one sigma error of < ±11.4° becomes fea- 
sible even with a modest acceptance of 0(10%). This er- 
ror can, of course, be reduced to the level of a few percent 
as the acceptance is improved. The method described to 
make use of inclusive states is likely to have wider appli- 
cation to the extraction of CP violating phases. 
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FIG. 1: The Xrnin vs - 8 f° r the toy model calculation given 
N B = 10 9 . The thin solid line is the result for D° -» K~-k + 
alone. The dashed line the result for CPES containing Ks 
together with related CPES containing Kl ■ The dotted lines 
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Case 


Accuracy 


CPES with K s and with K L 


±8.5° 


CPNES K-7T+ with K s and with K L 


±5° 


The CPNES K'n+ together with CPES, both 
with Ks only 


±9.0° 


K- + X together with K s CPES 


±2.5° 


K~~ +X together with K s as well as K L CPES 


±2.4° 



TABLE I: Attainable one sigma accuracy with various data 
sets given Nb = 10 9 ; note the 2nd and 5th cases are omitted 
from Fig 1 for clarity. 



Since overdetermining the system of equations is key 
to improving accuracy on 5 it may be useful to introduce 
additional constraints. First of all, one can replace the 
D° and the Ks with higher resonances which will tend 
to increase the total statistics. In addition, as suggested 
in ]l(| if the D° is replaced with a D ** then we can tag 
the flavor of the D** through the decay D ** — > D + n~. 
The analysis of decays with this tag thus reduces to that 
of There is the practical problem implementing 

this method that one must separate the D\ and Z?2 states 
which are 40 MeV apart. Secondly, as suggested in [fLs) , 
using the methods of (l^, [H]] one can directly determine 
Rf and f]£> from studies at a ^(3770) charm factory. 

Finally, the technique discussed here for replacing a 
single state with an inclusive one in the interference of 
two amplitudes has an immediate application to getting 
7. In |p| we consider this for the method of |l3| for 
extracting 7 from B~ — > K~D° with various D" final 
states. In particular, this gives a model independent way 
of analyzing three body B and D final states. 

In conclusion, we show that the ability to determine 8 
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